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Abstract 

The wobbling motion excited on triaxial superdeformed nuclei is studied in terms of the cranked 
shell model plus random phase approximation. Firstly, by calculating at a low rotational frequency 
the 7-dependence of the three moments of inertia associated with the wobbling motion, the mech- 
anism of the appearance of the wobbling motion in positive-7 nuclei is clarified theoretically — the 
^ ; rotational alignment of the ™ 13/2 q uasiparticle(s) is the essential condition. This indicates that the 

wobbling motion is a collective motion that is sensitive to the single-particle alignment. Secondly, 
we prove that the observed unexpected rotational-frequency dependence of the wobbling frequency 
is an outcome of the rotational-frequency dependent dynamical moments of inertia. 
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I. INTRODUCTION 



Deformation of the nuclear shape from spherical symmetric one has long been one of the 
most important issues in nuclear structure physics. Among them, searches for evidences of 
the triaxial (Y22 or 7) one have been pursued long time, for example, the even-odd energy 
staggering in the low-spin part of the 7 bands Q], the signature dependence of the energy 
spectra and the E2/M1 transition rates in medium-spin odd-odd and odd- A nuclei j^EHjL 
properties of the K isomers 0, |(| , and so on. But their results have not been conclusive; 
making a clear distinction between the static and the dynamic (vibrational) ones has not 
been successful up to now. Theoretically, appearance of the wobbling motion, which is 
well-known in classical mechanics of asymmetric tops j?| and whose quantum analog was 
discussed in terms of a rotor model about thirty years ago |8[, is a decisive evidence of 
static triaxial deformations. Subsequently its microscopic descriptions were developed by 
several authors 0, Q| • Since the small-amplitude wobbling mode carries the same quantum 
numbers, parity tt = + and signature a = 1, as the odd-spin members of the 7 band, Ref . \v\ 
anticipated that it would appear as a high-spin continuation of the 7 band, but it has not 
been resolved that in what nuclei, at what spins, and with what 7 wobbling modes would 
be observed. 

Shimizu and Matsuyanagi and Onishi performed extensive numerical calcula- 
tions for normally-deformed Er isotopes with relatively small |7|. Matsuzaki Q|, Shimizu 
and Matsuzaki and Horibata and Onishi also studied 182 0s with relatively large 
negative-7 but their correspondence to experimental information has not been very clear. 

These studies indicate the necessity of high-spin states in stably and strongly 7-deformed 
nuclei. Bengtsson studied high-spin states around 164 Hf [itJ and found systematic existence 
of the TSD (triaxial super- or strongly deformed) states with e 2 ~ 0.4 and I7I ~ 20°. This 
confirmed the discussion on the shell gap at N = 94 in Ref. ^| , the work in which the yrast 
TSD band in 163 Lu was reported; and in 2000 an excited TSD band was observed in this 
nucleus and from the strengths of the interband E2 transition rates this was unambiguously 
identified with the wobbling motion This data was analyzed by using a particle- rotor 
model [3| and the E2 transition rates were reproduced well. Subsequently TSD bands 
were found in some Lu and Hf isotopes and wobbling excitations were observed also in 



165,167 



Lu 



21 



A close look at these data, however, tells us that the sign of their 7- 



2 



deformation seems to contradict to an irrotational motion and that the unexpected behavior 
of the wobbling frequency has not been explained yet. 

Thus in the preceding Rapid Communication 23] we presented an answer to these prob- 
lems. In the present paper, after summarizing the discussion there we extend numerical 
analyses to elucidate it. An emphasis is put on the behavior of the calculated dynamic 
moments of inertia. 



II. WOBBLING MOTION IN TERMS OF THE RANDOM PHASE APPROXI- 
MATION 

We start from a one-body Hamiltonian in the rotating frame, 

ti = h- fuj TOt J x , (1) 

h = h m -A T (P} + P T )-\ T N T , (2) 
hmi = ^ + \M(u 2 x x 2 + u 2 y y 2 + u 2 z z 2 ) 

+Vls \-s + v n (\ 2 -(l 2 } No J. (3) 

In Eq.(J2J), t = 1 and 2 stand for neutron and proton, respectively, and chemical potentials A T 
are determined so as to give correct average particle numbers (N T ). The oscillator frequencies 
in Eq.(j3J) are expressed by the quadrupole deformation parameters e 2 and 7 in the usual way. 
They are treated as parameters as well as pairing gaps A T . The orbital angular momentum 
1 in Eq.Q is defined in the singly-stretched coordinates x' k = y^^fc, with k — 1 - 3 
denoting x - z, and the corresponding momenta. By diagonalizing h! at each cj rot , we obtain 
quasiparticle (QP) orbitals and the nuclear yrast (OQP) state. Since hi conserves parity n 
and signature a, nuclear states can be labeled by them. Nuclear states with QP excitations 
are obtained by exchanging the QP energy and wave functions such as 

(-e^V^Lg^^U^Vp), (4) 

where p, denotes the signature partner of ji. 

We perform the random phase approximation (RPA) to the residual pairing plus doubly- 
stretched quadrupole-quadrupole (Q" ■ Q") interaction between QPs. Since we are interested 
in the wobbling motion that has a definite quantum number, a = 1, only two components 
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out of five of the Q" ■ Q" interaction are relevant. They are given by 



(-) 

hit 



(5) 



if =1,2 



where the doubly-stretched quadrupole operators are defined by 

Q'k = Qi<{ x k — > x'l = — X k), 



and those with good signature are 



1 



(Q K ± Q- K ) . 



(6) 



(7) 



%1 + Sko) 

The residual pairing interaction does not contribute because P T is an operator with a = 0. 
The equation of motion, 



h' + H<-\xt 



J RPA 



fkj n Xl, 



for the eigenmode 



(a=±l/2) 



il<V 



leads to a pair of coupled equations for the transition amplitudes 



T, 



K,n 



Q ( k\K 



(8) 



(9) 



(10) 







Then, by assuming 7 ^ 0, this can be cast [lCj into the form 



Jx~Jy \^n) ) ( Jx — Jz 



(off). 



0. 



(11) 



which is independent of k> k ^s. This expression proves that the spurious (Nambu-Goldstone) 
mode given by the first factor and all nomal modes given by the second are decoupled from 



each other. Here 



in Refs 



.0 



(J x )/uj TOt as usual and the detailed expressions of Jif (u n ) are given 



6 \J x 

l4 Il5l | . Among normal modes, one obtains 

{^Jx — Jy ^(^wob)^ {^Jx — Jz ^(^wob) 



wob k'rot 



(12) 



Jy \^ob)<Jz ^(u; wo b) 

by putting uj n = c<j wo b. Note that this gives a real excitation only when the right-hand side 
is positive and it is non-trivial whether a collective solution appears or not. Evidently this 



coincides with the form derived by Bohr and Mottelson in a rotor model 8] and known 

Q 

in classical mechanics |7], aside from the crucial feature that the moments of inertia are 
cjrot-dependent in the present case. 

One drawback in our formulation is that our J x tends to be larger than corresponding 
experimental values because of the spurious velocity dependence of the Nilsson potential 



as discussed in Refs. 



24, 25]. as discussed there but that for Jy^ 1 has 

not been devised yet. Therefore we assume for the present a similar discussion holds for 
the latter and accordingly the ratio Jy*P (uJ wo b) / J x which actually determines oo woh is more 
reliable than their absolute magnitudes. 

Interband electric quadrupole transitions between the n-th excited band and the yrast 
are given as 



in terms of 



B(E2 : I B - (J ± l) yrast ) = \ {Ti E J ± T^)\ (13) 



= ejT K , n . (14) 



They will be abbreviated to B(E2) ont later for simplicity. In-band ones are given as 

B(E2 : 1 I - 2) = \ (Q<+>< £ >) + \ {Q'^')] , (15) 



in terms of 



of"") =4 <«?')■ < ie > 



and assumed to be common to all bands. They will be abbreviated to B(E2) in . Here we 
adopted a high-spin approximation 2f|. The transition quadrupole moment Q t is extracted 
from B(E2) in by the usual rotor-model prescription. 

To compare collectivities of these two types of E2 transitions, we introduce a pair of 
deformation parameters 

Then it is evident that the in-band one is expressed as 




B(E2 : / -> / - 2) = l -R' - a ™)\ (18) 



As for the interband ones, by expanding ^ by X^s and X n s, where n runs both normal 



modes and the Nambu-Goldstone mode X^ G = -^j(J z + iJ y ), we obtain from [Q[ \q[ ^ 
a kind of sum rule 



^ Tl, n T 2) n = - J R 4 OL y OL z 



(19) 



n^NG 



Consecutively introducing the ratios of the dynamic to static deformations, 



y,n 



2R 2 



T2,n 



2R 2 



the sum rule above reads 



r y,n r z,n 



(20) 



(21) 



n^NG 



The dynamic amplitudes Tk, u describe shape fluctuations associated with the vibrational 
motion in the uniformly rotating frame. Transformation to the body-fixed (Principal-Axis) 
frame turns the shape fluctuation into the fluctuation of the angular momentum vector, 
i.e., the wobbling motion. This transformation relates the ratios, r y>n and r ZjU , to the 
moments of inertia |15j |: 



1 



w g , n 



21 \W, 



1/4 



yji 



1/4 



— , , (22) 

where c n is a real amplitude that relates the dynamic amplitude Tif, n and the moment of 
inertia, a n is the sign of \ J X — Jy^\ (so o n > for wobbling-like RPA solutions), and 



W z>n = l/jW(u n )-\/J x . 



(23) 



Thus, the interband B(E2) is rewritten as, 



B(E2 : I n -> (7±r 



yrast J 



4 „2 



a 



1/1 



IV \ 1/4 ' 



1 2 



(24) 
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which coincides with the formula given by the rotor model |8(, except for the appearance of 
the amplitude c n and sign o n . Substituting the ratios, r y ^ n and r Zj7l , into Eq. (}2"T]) . one finds 
that the amplitudes should satisfy 

n^NG 

This form of sum rule clearly indicates that the amplitude c n is a microscopic correction 
factor quantifying the collectivity of the wobbling motion, for which c\ ~ 1 means the full- 
collectivity and reproduces the results of the macroscopic rotor model in both the energy 
and the interband B(E2) values. 

III. NUMERICAL CALCULATION AND DISCUSSION 
A. Summary of the preceding study 

Since the first experimental confirmation of the wobbling excitation in 163 Lu 7 ~ 
+20° has been widely accepted as the shape of the TSD states in this region. This is 
predominantly because the calculated energy minimum for 7 ~ +20° is deeper than that 



for 7 ~ —20° |17j according to the shape driving effect of the aligned 7rii 3 / 2 quasiparticle. 
The recent precise measurements of Q t [27] also support this. On the other hand, the 
sign of 7-deformation leads different consequences on moments of inertia, which are directly 
connected to the excitation energy of the wobbling mode through the wobbling frequency 
formula p], c.f. Eq. (jl2J) . Since the RPA is a microscopic formalism, no distinction between 
the collective rotation and the single-particle degrees of freedom has been made. 

Therefore, the moments of inertia calculated in our RPA formalism in sect]n]are those 
for rotational motions of the whole system. In contrast, the macroscopic irrotational-like 
moments of inertia are often used in the particle-rotor calculations, where J y > J x ^> J z 
for 7 ~ +20° and they lead to an imaginary wobbling frequency co> wo b- It is, however, noted 
that the moments of inertia of the particle-rotor model are those of the rotor and no effect 
of the single-particle alignments is included, so that they do not necessarily correspond to 
those calculated in our RPA formalism. 

In the preceding paper j2j] we have performed microscopic RPA calculations without 
dividing the system artificially into the rotor and particles. That work proved that for the 
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calculated moment of inertia, J x = (J x )/u} TOt , the contribution from the aligned QP(s), 
AJ X = iQp/ti; r ot with zqp being the aligned angular momentum, is superimposed on an 
irrotational-like moment of inertia (J~ y > <J X ) of the "core". Consequently the total J x is 
larger than J y , which makes wobbling excitation in 7 > nuclei possible. 

The second consequence of the formulation adopted in Ref. [23( is that the three moments 
of inertia are automatically uj rot - dependent even when the mean-field parameters are fixed 
constant. This is essential in order to explain the observed u; ro t-dependence of u; wo b - 
decreasing as co> rot increases. Otherwise o; wob is proportional to c<j rot . 

Another important feature of the data is that the interband B(E2) values between the 
wobbling and the yrast TSD bands are surprisingly large. Our RPA wave function gave 
extremely collective B(E2) out that gathered |c n=wo b| — 0.6 - 0.8 in the sum rule (Eq. (|25j) ) 
but the result accounted for only about one half of the measured one. 

To elucidate these findings more, in the following we extend our numerical analyses 
putting a special emphasis on the 7-dependence of the moments of inertia in subsect flll Bl 
Dependence on other parameters is also studied in detail. Features in common and different 
between even-even and odd- A nuclei are also pointed out. In subsect flll C\ we discuss uo rot - 
dependence. In subsect lTlI Dl characteristics of B(E2) out are discussed. Calculations are 
performed in five major shells; N osc = 3-7 for neutrons and N osc = 2-6 for protons. The 
strengths Vi s and Vu in Eq.(j2J are taken from Ref. . 



B. Dependence on the mean-field parameters 7, €2, and A 

1. The even-even nucleus ie8 Hf 

n 

Hafnium-168 is the first even-even nucleus in which TSD bands were observed [23]. In 
this nucleus three TSD bands were observed but interband 7-rays connecting them have 
not been observed yet. This means that the character of the excited bands has not been 
established, although we expect at least one of them is wobbling excitation. An important 
feature of the data is that the average transition quadrupole moment was determined as 
Qt = 11.4^2 eb. This imposes a moderate constraint on the shape. Referring to the weak 
parameter dependence discussed later, we choose e 2 = 0.43, 7 = 20°, and A n = A p = 0.3 
MeV, which reproduce the observed Q t , as a typical mean-field parameter set. 



8 



First we study the dependence of various quantities on 7 and other mean-field parameters 
at hu> rot = 0.25 MeV. Around this frequency the (vrz 13 / 2 ) 2 alignment that is essential for 
making wobbling excitation in 7 > nuclei possible is completed and therefore the wobbling 
motion is expected to emerge above this frequency (see FigJ7| shown later). 

Figure H] shows dependence on 7 calculated with keeping e 2 = 0.43 and A n = A p = 0.3 
MeV. Figure [T^a) graphs the calculated excitation energy in the rotating frame, fta> W ob- As 
7 comes close to (symmetric about the z axis) and —60° (symmetric about the y axis), 
^wob approaches 0, see Eq.(fT2"j). We did not obtain any low- lying RPA solutions at around 
7 = 40° whereas a collective solution appears again for 50° < 7 < 60°. 
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FIG. 1: Triaxiality dependence of (a) excitation energy of the wobbling motion, (b) three moments 
of inertia associated with it, (c) expectation values of angular momenta in the yrast state, and (d) 
quadrupole transition amplitudes between the wobbling and the yrast states in 168 Hf, calculated 



at huj 



rot 



0.25 MeV with e 2 = 0.43 and A n = A„ = 0.3 MeV. 



Figure ^b) shows the calculated moments of inertia. Their 7-dependence resembles the 
irrotational, the so-called 7- reversed, and the rigid-body moments of inertia, in 7 < 0, 
< 7 < 40°, and 50° < 7 < 60°, respectively. These model moments of inertia are given by 
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Jl" = 4Bf3 2 sin 2 (j + ^irk), 



(26) 



and 



JT = W 2 sin 2 (- 7 + ^A;), (27) 



^ riS = ^o(l-J^cos( 7 + ^)], (28) 



where k = 1 - 3 denote the x - z principal axes, B the irrotational mass parameter, J§ the 

rigid moment of inertia in the spherical limit, and (3 is a deformation parameter like 62. The 

7-reversed moment of inertia was introduced to describe positive-7 rotations in the particle- 
fl 

rotor model [3J but its physical meaning has not been very clear; in particular, it does 
not fulfill the quantum-mechanical requirement that the rotations about the symmetry axis 
should be forbidden. We have clarified in the preceding paper [23] that the contributions from 
aligned quasiparticles superimposed on irrotational-like moments of inertia (J x < J y ) can 
realize J x > J y and this is the very reason why the wobbling excitation (see Eq.(|12Jl) appears 
in positive-7 nuclei. We also discussed that multiple alignments could eventually lead to a 
rigid-body-like moment of inertia. Figure ^c) indicates that, in the present calculation in 
which configuration is specified as the adiabatic quasiparticle vacuum at each co> rot , two 7Ui 3 / 2 
protons align for 7 > as mentioned above while they have not fully aligned for 7 < at this 
u rot . In other words, these figures cover the both regions in which the (vrz 13 /2) 2 alignment 
is necessary (7 > 0) and that is not necessary (7 < 0) for obtaining wobbling excitations. 
This aligned angular momentum determines the overall 7-dependence of J x in Fig^b). As 
for the neutron part, corresponding to the disappearance of the solution at around 7 = 40°, 
the expectation value of the neutron angular momentum, (J x ) n , drops around this region. 

To look at this more closely, we investigate the Nilsson single-particle diagram at u rot = 0. 
Figure Efa) graphs neutron single-particle energies for < e 2 < 0.43 with 7 = 0, while 
FigEfb) for < 7 < 60° with e 2 = 0.43. The chemical potential that gives correct neutron 
number N = 96 for 7 > at huj vot = 0.25 MeV is also drawn in the latter. This figure 
clearly shows that with this e 2 a shell gap exists for 7 < 20° at iV = 96. And by comparing 
this with FigHJwe see that the dropping of (J x ) n is a consequence of the deoccupation of the 
orbital that is [651 1/2] at 7 = (hereafter simply referred to as the [651 1/2] orbital even at 
7 7^ 0) originating from the mixed (#9/2-211/2) spherical shell. Figure Efb) also explains the 

10 



reason why the wobbling excitation revives at around 7 = 50° again; the occupation of other 
oblate-favoring orbitals such as [503 7/2] makes it possible and leads to a rigid-body-like 
behavior of the moments of inertia. Figures |5Jc) and (d) are corresponding ones for protons. 
This indicates that the proton shell gap is robuster. 
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FIG. 2: Nilsson single-particle energy diag ranis cit u ro t — 

0, (a) for < e 2 < 0.43 with 7 = 0, and 
(b) for < 7 < 60° with €2 = 0.43 for neutrons, (c) and (d) are corresponding ones for protons. 
Solid and dashed curves represent even and odd parity orbitals, respectively. Asymptotic quantum 
numbers of some important orbitals are explicitly indicated. Chemical potentials that give particle 
numbers N = 96 and Z = 72 for 7 = +20° at Hu ro t = 0.25 MeV are also indicated in (b) and (d). 

Figure Hfd) graphs the quadrupole transition amplitudes Tk (K = 1,2) associated with 
the wobbling mode. (T K corresponds to (— 1) k ~ 1 Qk in RefJl5|L_) This shows that their 

151 ] . This feature can be 



14, 



relative sign changes with that of 7 as discussed in Refs. 
understood as: 7 ~ is the 7-vibrational region because the K = 2 component is dominant 
(see also J x ~ Jy and — in FigJTfb)), and the mixing of the K — 1 component due 
to triaxiality and rotation gives rise to the character of the wobbling motion. This relative 
sign leads to a selection rule of the interband transition probabilities B(E2) ont | -3J . In the 
present case we obtain B(E2 : I —> I — l) out ^ B(E2 : 1 — > 1+ l) out for 7^0, and typically 



11 



their ratio to the in-band ones is B(E2 :/—>•/— l) oat /B(E2 : I -> I - 2) in ~ 0.1. 

Figure El shows dependence on e 2 calculated with keeping 7 = 20° and A n = A p = 0.3 
MeV. The steep rises at around e 2 = 0.33 in FigsJ^a) and (b) indicate the necessity of the 
(7rz 13 / 2 ) 2 (the [660 1/2] orbital in FigO^c)) alignment for the appearance of the wobbling 
mode although the critical value of e 2 itself is frequency dependent. Aside from this, u; wo b is 
almost constant in the calculated range. The slight increase at around e 2 = 0.4 stems from 
the occupation of the z/[651 1/2] orbital. We have confirmed that in this case the (^jis/ 2 ) 2 
alignment at around e 2 = 0.47 seen in Fig|3fb) does not affect u; wob visibly since AJy eS ^ in 
this case is almost the same as AJ X although its reason is not clear. Figure Etc) graphs Q t - 
This figure indicates that the chosen shape e 2 = 0.43 and 7 = 20° reproduces the measured 
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FIG. 3: Deformation dependence of (a) excitation energy of the wobbling motion, (b) expectation 
values of angular momenta in the yrast state, and (c) transition quadrupole moment in the yrast 
state in 168 Hf, calculated at huj TOt = 0.25 MeV with 7 = 20° and A n = A p = 0.3 MeV. 



Figure Ufa) shows dependence on the pairing gaps. Since we do not have detailed in- 
formation about the gaps, we assume A n = A p for simplicity. This figure shows that the 
dependence on the gaps is weak unless they are too large. Since the static pairing gap A 
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is expected to be small, say, A < 0.6 MeV, in the observed frequency range, c<j wo b is not 
sensitive to the value of A. This is a striking contrast to the (3 and 7 vibrations; it is well 
known that pairing gaps are indispensable for them. Here we note that the behavior of the 
^wob correlates well with (J x ) p presented in Fig0Jb). 
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FIG. 4: Pairing gap dependence of (a) excitation energy of the wobbling motion and (b) expectation 
values of angular momenta in the yrast state in 168 Hf, calculated at huj ro t = 0.25 MeV with e<i = 0.43 
and 7 = 20° . A n = A p is assumed for simplicity. 



2. The odd- A nucleus e 'Lu 

Next we study 167 Lu in a way similar to the preceding 168 Hf case. We choose 7 = 20° and 
A n = A p = 0.3 MeV as representative mean-field parameters as above. As for €2, however, 
we examined various possibilities because Q t has not been measured in this nucleus. Since 
the sensitive e2-dependence through the occupation of the [651 1/2] orbital appears only at 
fy^rot > 0.4 MeV and therefore the "band-head" properties do not depend on €2 qualitatively, 
first we discuss them adopting €2 = 0.43 in order to look at the difference between the even- 
even and the odd-Z cases. 

Figure El shows dependence on 7 at hu ro t = 0.25 MeV with keeping e 2 = 0.43 and 
A n = A p = 0.3 MeV constant. Figure a) graphs tu wo b- In the 7 > region, the solution 
is quite similar to the 168 Hf case. In the 7 < region, that for —60° < 7 < —30° is quite 
similar again but for —30° < 7 < its character is completely different. In this region the 
presented solution is the lowest in energy and becoming collective gradually as 7 decreases. 
The largeness of cj wo b corresponds to that of J x — Jy CS ^ in FigEl^b). Comparison of FigsEf c) 
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and njc) certifies that the alignment of the irii3/2 quasiparticle(s) is almost complete for 
7 > whereas less for 7 < 0. This produces quantitative even-odd differences as explained 
below. 




Having confirmed that these features are independent of 62 and N except that we did 
not obtain any low-lying solutions for 35° < 7 < 60° in the small-e2 cases, we look into 
underlying unperturbed 2QP energies to see the even-odd difference. In FigEl we present 
the energies of the lowest (7r(iV osc = 6)) 2 states which represent the biggest difference. In 
the yrast {niiz/2) 2 configuration, A p and B p in the usual notation are occupied in the even-Z 
case, the lowest 2QP state of signature a = 1 with respect to this is B P C P (where ~~ denotes 
the conjugate state, see Eq.(0}). In the odd-Z case in which A p is occupied, the lowest one is 
B p Ap. Since both e' Bp and decrease as 7 decreases, this 2QP state becomes the dominant 
component in the lowest-energy RPA solution. Note here that the sum e' Bp corresponds 
to the signature splitting between A p and B p when they are seen from the usual even-even 
vacuum. Since both B p and A p are of K = 1/2 character, the resulting RPA solution can 
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not have the K = 2 collectivity as shown in Figl^d). According to the relation 



115], 



Jt 



(eff) 



1 + 



^wob 



SU17 



uj mt sin (7 + hz) T 2 



(29) 



J, (eff) in Fig^b) becomes small for —30° < 7 < 0. These discussion serves to exclude the 
possibility of 7 ~ —20° for the TSDs that support collective wobbling excitations in the 
odd-Z cases, whereas the even-odd difference in 7 > is merely quantitative. 
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FIG. 6: Energies of the lowest (7r(N osc = 6)) 2 two quasiparticle states in 168 Hf and 1B7 Lu, calculated 



167t 



at the same time as FigsQandEl respectively. 



C. Dependence on the rotational frequency w ro t 
1. 168 Hf and 174 Hf 

The analyses above indicate that the chosen mean-field parameters are reasonable, and 
therefore we proceed to study cj rot -dependence with keeping these parameters constant. 
Figure E| shows the result for 168 Hf. These figures indicate again the (vri 13 / 2 ) 2 alignment that 
makes J^ x larger than J~y eS ^ is indispensable for the formation of the wobbling excitation. At 
around hw rot = 0.45 MeV the (VJ15/2) 2 alignment occurs. In contrast to the low-frequency 
case reported in FigJSl in the present case its effect on cu wo b is visible as a small bump. 



Dands has not been resolved, some 
We suggest this is related to the 



Although the character of the observed excited TSD 
anomaly is seen at around this u rot in one of them 2 
(z/j 15 / 2 ) 2 alignment since this is the only alignable orbital in this frequency region of this 
shape. However we note that in 167 Lu an interaction with a normal deformed state at around 



this frequency is discussed in Ref. 



22|- 
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We performed calculations also for 7 = —20°. In that case, however, wobbling excitation 
exists only at small u ro t because J x — is small as seen from FigHJ^b). 
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FIG. 7: Rotational frequency dependence of (a) excitation energy of the wobbling motion, (b) ex- 
pectation values of angular momenta in the yrast state, and (c) three moments of inertia associated 
with the wobbling motion in 168 Hf, calculated with e 2 = 0.43, 7 = 20° and A n = A p = 0.3 MeV. 



Very recently TSD bands were observed in another even-even nucleus, 174 Hf <:>U|. it is 
not trivial if a similar band structure is observed in the nucleus with six neutrons more 
since the existence of the TSD states depends on the shell gap. Multiple TSD bands were 
observed but connecting 7-rays have not been resolved also in this nucleus. We performed a 
calculation adopting e 2 = 0.453 and 7 = 16° suggested in Ref. Q and A n = A v = 0.3 MeV. 
The result is presented in Fig|H] The most striking difference from the case of 168 Hf above 
is that u; wob decreases steadily as u rot increases after the (vrii 3 / 2 ) 2 alignment is completed. 
This is because the (z/j'15/2) 2 alignment that causes the small bump in the 168 Hf case shifts 
to very low uo rot due to the larger neutron number. 
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FIG. 8: Rotational frequency dependence of excitation energy of the wobbling motion in 174 Hf, 
calculated with e 2 = 0.453, 7 = 16° and A n = A p = 0.3 MeV. 

2. W7 Lu 

The wobbling excitation was first observed experimentally in 163 Lu [3], later it was also 
observed in 165 Lu jsjj and 167 Lu (2^]. The characteristic features common to these isotopes 
are 1) u; wo b decreases as u TOt increases contrary to the consequence of calculations adopting 
constant moments of inertia, and 2) B(E2 — l) out /B(E2 — 2) in is large - 

typically around 0.2. 

Here we concentrate on the isotone of 168 Hf discussed above, that is, 167 Lu in order to see 
the even-odd difference. A comparison of Figs|7| and El proves that all the differences are due 
to the fact that the number of the aligned 7rii3/2 quasiparticle is less by one: 1) The 
alignment at around hjj TO t = 0.2 MeV is absent, and 2) the B P C P crossing occurs at around 
hu> rot = 0.55 MeV, which is proper to the configuration. Figure 0a) shows that 

our calculation does not reproduce the data, although in each frequency range in which the 
configuration is the same u; wo b decreases at high u rot as in the cases of the even-even nuclei 
presented above. This result might indicate that there is room for improving the mean field. 
The J x in Fig|2fc) is larger than the experimentally deduced value by about 20 - 30 %. 
This is due to the spurious velocity dependence of the Nilsson potential mentioned in SeclTTl 



D. Interband B(E2) transitions 

Compared to the excitation energy, the interband B(E2) values relative to the in-band 
ones have been measured in only few cases. In FiglT0| we report calculated B(E2) ratios 
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FIG. 9: The same as Fig|7|but for 167 Lu. Experimental values taken from Ref. 22] are also included 
in (a). 

for / (wobbling on yrast TSD) — > I ± 1 (yrast TSD) transitions in 168 Hf and 167 Lu. The 
measured ones are also included for the latter. 

The first point is the magnitude of the larger (/ — > 1 — 1) ones. Apparently, the calculated 
B(E2) values are smaller by factor 2-3. The measured interband B(E2) values amount 
almost to the macroscopic rotor value. In the RPA calculations, as summarized in sect 01 the 
B(E2) value is reduced by a factor c^ =wob (see Eq. fj2lj)): Only the case with the full-strength 



-'n=wob 



1 the rotor value is recovered. Although the obtained RPA wobbling solutions are 



extremely collective in comparison with the usual low-lying collective vibrations, like the 
(3- or 7-vibrations, for which typically \c n \ ~ 0.3 - 0.4, this factor is still |c„= wo b| — 0.6 
- 0.8. This is the main reason why the calculated B(E2) values are a factor 2 - 3 off 
the measured ones. As is well-known, giant resonances also carry considerable amount of 
quadrupole strengths, so it seems difficult for the microscopic correction factor c^ =wob to be 
unity; it is not impossible, however, because the "sum rule" discussed in sectOJis not the 
sum of positive-definite terms. In the RPA formalism, the reduction factor c^ =wob for the 
B(E2) value, Eq. (J2IJ), comes from the fact that the wobbling motion is composed of the 
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coherent motion of two-quasiparticles, and reflects the microscopic structure of collective 
RPA solutions. The measurement that the B(E2) value suffers almost no reduction may be 
a challenge to the microscopic RPA theory in the case of the wobbling motion. Calculated 
B(E2) ratios for 174 Hf are slightly smaller than those for 168 Hf in Fig lTUT a). 
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FIG. 10: Interband E2 transition rates for / (wobbling on yrast TSD) — ► I ± 1 (yrast TSD) 
transitions in (a) 168 Hf and (b) 167 Lu. The latter is presented as functions of 2x spin /, while 
the former is presented as functions of the rotational frequency since experimental spin assignment 
has not been done for 168 Hf. The rotational-frequency range corresponding to (b) is very narrow 
in comparison to (a). Interband transition rates are divided by the in-band ones. Experimental 
values zj] are also shown in (b). Noting that, for 167 Lu, the states I + 1 (TSD1) are slightly higher 
in energy than I (TSD2) at I > 5l/2h and B(T X ;I -»■ I + 1) ~ B(T X ;I + 1 -»■ I) holds at high 
spins, we plotted those for / — > / + 1 at the places with the abscissae I + 1 in order to show clearly 
their characteristic staggering behavior. 



The second point is the staggering, that is, the difference between / — > / ± 1. We 
clarified |3| its unique correspondence to the sign of 7 as mentioned in subsect llHBI that 
holds for both even-even and odd-A systems. Recently this staggering was discussed from a 
different point of view jsi|; but it looks to apply only to 7 < cases. 



IV. CONCLUSION 



The nuclear wobbling motion, which is a firm evidence of stable triaxial deformations, 
was identified experimentally in the triaxial superdeformed odd- A Lu isotopes. In principle, 
wobbling excitation is possible both in 7 > and 7 < nuclei. Every information, theoretical 
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and experimental, suggests 7 > for these bands. According to the wobbling frequency 
tormui : g ,, E, ra , H, e XClt a tl on in ouclei rotate Pri n clp ; ily aW th e a; axis requires 
Jx > JyiJz-, although irrotational-like model moments of inertia give J x < J y for 7 > 
0. To solve this puzzle, we studied the nuclear wobbling motion, in particular, the three 
moments of inertia associated with it in terms of the cranked shell model plus random 
phase approximation. This makes it possible to calculate the moments of inertia of the 
whole system including the effect of aligned quasiparticle(s). The results indicate that the 
7-dependence of the calculated moment of inertia is basically irrotational-like (3 X ^ J y for 
7 5s 0) if aligned quasiparticle(s) (vrii 3 / 2 in the present case) does not exist. But once it is 
excited, it produces an additional contribution, AJ X = iQp/ui rot , and consequently can lead 
to J x > J y . This is the very reason why wobbling excitation exists in 7 > nuclei. In 
this sense, the wobbling motion is a collective motion that is sensitive to the single-particle 
alignments. 

The resulting moment of inertia for < 7 < 30° resembles the 7- reversed one, i.e., the 
irrotational moment of inertia but with J x and J y being interchanged. That for 50° < 7 < 
60°, where single-particle angular momenta dominate, is rigid-body- like. That for 7 < 
is irrotational-like except for odd-v4 nuclei with —30° < 7 < where a specific 2QP state 
determines the lowest RPA solution. 

Having studied qualitative features of the three moments of inertia at a low rotational 
frequency, we calculated wobbling bands up to high uj TOt . Experimentally they were con- 
firmed only in odd- A Lu isotopes as mentioned above. The most characteristic feature of the 
data is that u; wo b decreases as u) rot increases. This obviously excludes constant moments of 
inertia. In our calculation three moments of inertia are automatically to> ro t-dependent even 
when mean field parameters are fixed constant. It should be stressed that the wobbling-like 
solution in our RPA calculations is insensitive to the mean-field parameters, especially to 
the pairing gaps, as is shown in subsect lill B 11 This distinguishes the wobbling-like so- 
lution from the usual collective vibrations, which are sensitive to the pairing correlations. 
Thus, our microscopic RPA calculation confirms that the observed band is associated with 
a new type of collective excitation, although comparisons to the observed excitation energy 
indicate that there is room for improving the calculation. 

As for the interband transition rates, our calculation accounted for only about one half 
or less of the measured ones, even though the wobbling-like solution is extremely collective 
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compared to the usual vibrational modes. This issue is independent of the details of choosing 
parameters. This confronts microscopic theories with a big challenge. 
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